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The interaction between gravitational and electromagnetic waves in the presence ol a static mag- 
netic field is studied. The field strength of the static field is allowed to surpass the Schwinger critical 
field, such that the quantum electrodynamical (QED) effects of vacuum polarization and magneti- 
zation are significant. Equations governing the interaction are derived and analyzed. It turns out 

f~*) ' that the energy conversion from gravitational to electromagnetic waves can be significantly altered 

7-H ' due to the QED effects. The consequences of our results are discussed. 

o: 

^ . PACS numbers: 04.30.Nk, 12.20.Ds, 04.30.Tv 

>>: 

1^ . I. INTRODUCTION 

OO ■ As studied by many authors jl'-l.l6| there exist numerous mechanisms for the conversion between gravitational waves 

Cn (GWs) and electromagnetic (EM) waves. In particular, the propagation of GWs across an external static magnetic 

■ field gives rise to a linear coupling to the electromagnetic field (see e.g. Ref. [ll-Q), which may lead to the GW 

O ] excitation of ordinary EM waves in vacuum, or of roagnetohydrodynamic (MHD) waves in a plasma [J-Q. Many 

^ • nonlinear coupling mechanisms are also possible [J, I6l4l0l| . Conversion of energy from gravitational to electromagentic 

5^ ' degrees of freedom has been pointed out as a means to indirect detection of gravitational waves by several authors 

{^, (see e^. Refs. SSllll), since the latter is so much easier to detect. For astrophysical application (see e.g. Refs. 

Ul, B M, M \1M, llfll)) naturally this requires well developed theories to recognize the signature of the gravitational 

I , origin. Furthermore, there must be a sufficient amount of energy conversion taking place. Specifically, considering 

^ ■ the coupling due to a static magnetic field, it has been noted that more energy can be converted from gravitational to 

^sD i electromagnetic degrees of freedom if the interaction region is larger, and if the magnitude Bq of the static magnetic 

OO '. field is larger 3] . In the case that the interaction region is magnetized vacuum, with a size smaller than the background 

^^ curvature radius, it has been found that the energy converted is linear in the background field energy de nsity ^2, 3]. 

This result, however, does not account for quantum electrodynamic (QED) vacuum polarization effects |17l - [20| . which 

lO ' become significant when Bq approaches the value Ecr/c, where Ecr = mlc^ /he ~ 10"'^*V/m is the Schwinger critical 

^3^ , field, TTie is the electron mass, e is the elementary charge, c is the speed of light in vacuum, and h = 2Trh is the Planck 

constant. 

In the present paper we will investigate the QED influence on gravitational-electromagnetic interaction in a static 
field Bq that may be stronger than the characteristic QED scale Ea/c. It should be noted that such intense field 
do occur in nature, specifically close to magnetars where close to the surface the magnetic field strength may reach 
10^*^ — lO^-'^T f2l|. Starting from Einstein's equations, together with the Heisenberg-Euler Lagrangian to describe 
C^ ' vacuum polarization and magnetization in the electromagnetic theory, the basic equations for small amplitude wave 
propagation on a background with a strong static magnetic field Bq is derived. In order to simplify the calculation, 
the size of the interaction region is assumed to be much smaller than the background curvature. It is found that 
the vacuum polarization effects lead to a saturation, such that the energy conversion (almost) stops to grow with Bq 
beyond a certain value -Bsat- This value depend on the length L of the interaction region. For a large L, the saturation 
value is much smaller than the QED scale, i.e. i?sat ^ -E-cr/c (in which case the weak field QED corrections [13| of 
the Heisenberg-Euler theory would have sufficed), but for shorter interaction regions we may have E„/c <C Ssat in 
which case the full theory is required. The relevance of our model calculation to astrophysical problems is discussed 
at the end of the paper. 
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II. BASIC EQUATIONS 

The Lagrangian for soft photon (i.e. photon energy much smaller than electron rest mass energy) light propagation, 
taking one loop corrections into account, is given by |19l. [20l. [23| 



1 T- " r°° ^^ -eE s/c 

Mo 27rpioe^ Jo « 



(es)^a6 coth [eas) cot [ehs) - (a^ - b^) 



Ao.r , (1) 



where a = [^(T^T^ + -F]i/2, h = [^^jT^^ _ j-ji/2 ^ ^ ^ (l/4)^„0i^"'3 , g ^ (1/4)F„^F"^ , F'^^ is the 
electromagnetic field tensor, F"'^ = ga/SM^^ii^ ^q^a''' ^j-^g totally antisymmetric tensor, Aq the four-potential, j" 
the four-current and a the fine structure constant. The Euler-Lagrange equations of motion for the Lagrangian ([1]) 
becomes 

T^^r + IqF^ + \ b^^F^^'F^P + 7ee^^'^^o/5]i^,f + 7^e [i^'^'^i^a/^ + ^''"^a/?] i^.^l^ = -f , (2) 

where we have applied the Eq. (2) of Ref. [1^ to a curved background, and introduced the quantities 

_ ac _dC 
^^^ = Wdg ■ ^^^ 

The physics of strong field vacuum polarization and vacuum magnetization is thus encoded in the parameters intro- 
duced in Eq. ([3]). For the case of interest to us, i.e. no external electric field, the scalars,7jF, 7g,7j^.F,7a0 ^^'^ Iro 
can be computed analytically as functions of the external constant magnetic field strength Bq. This procedure which 
involves the solution of numerous integrals is described in Ref. 'l9| , and the explicit expressions of the scalars can be 
found in Appendix A. 

In the paper we will study the infiuence of a GW on a strong magnetic field. The metric of a linearized GW 
propagating in the z-direction can be written 

ds^ = -c^df + (1 + h+) dx^ + (1 - h+) dy^ + 2hy,dxdy + dz^ (4) 

where the two independent polarizations /i+ and ft-x depend on the coordinates as /i+.x = /j+,x {z — ct). Furthermore, 
we define an orthonormal tetrad by 



eo = 


c 






ei = 


/ 1 ^ 


a,- 


- -^hy^dy, 


62 = 


^ (-K 


]dy- 


~\h.d., 



e3 - a, . (5) 

In linearized theory of gravity, the relevant components of the Einstein equations read: 

{el - dl)h+ = K(,5rn - 5T22). {el - dl)h^ = 2n{5Ti2), (6) 

where k = SttG/c , and G is the gravitational constant. The energy- momentum tensor associated with the Lagrangian 
^ is written T^,^ — —jj^F"Fai, + {Q^g — C^^Qiiv, see [23|, and expressions for (5Tii, 8T22 and ^Tyi, linearized around 
the strong magnetic field Bq, is worked out in Appendix A. 

Next we follow the covariant approach presented in Ref. [24] for splitting the EM and material fields in a 1 -I- 3 
fashion. Suppose an observer moves with 4- velocity m". This observer will measure the electric and magnetic fields 
Fa = FajiU^ and Ba = e^p^F^^ /2^ respectively, where Fap is the EM field tensor and Cq.^^ is the volume element 
on hyper-surfaces orthogonal to m". We also define the spatial gradient operator as V — (ei, 62, 63). Using the 1 -t- 3 



split we write the Maxwell equations in the tetrad basis ([5|). From Eq. ([2]) and the Faraday equation, F^ij.^^ = 0, we 
obtain 

(7) 

(8) 

(9) 

Je] , (10) 
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where jg is the combined vacuum polarization and vacuum magnetization current density, which from Eq. ^^ can 
be seen to take the form 

3q = -77^ ('2££f'='f'- + H^f'^'fA e,Fuu (11) 

^ 2/io V 7F IF J 

and the effective (i.e. gravity induced) charge densities and current densities are 

PE ^ -eo [ll^E^ + e^^-'llpcB^] , 
PB ^ -eo [7|.cB^ - t'^^''l%E,] , 

-hop - l$o)— + lip— - ^"^'h'opB, + I'p.Bs) 
-(7o/3 - 7/30)-^ + 70/3-^ + e '^ ' y°f^~ ^f^^~ ) 
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where the Greek indices takes values between 1 and 3, and the Latin indices between and 3. From here on we will 
be concerned with a GW wave propagating across a magnetic field. Explicit expressions of the source terms for this 
case is obtained by substituting the QED-parameters from Appendix A into Eq. (|TT]) . and the rotation coefficients 
for a linearized GW presented in Appendix B into Eq. (fT2|) . 



III. WAVE INTERACTION 

The most efficient interaction of a GW with a static magnetic field occurs if the GW propagates perpendicular 
to the magnetic field. As has been found by e.g. Refs. [2, 0], the fact that the GW fulfills the same dispersion 
relation as EM- waves, makes the energy conversion resonant. As a consequence, the energy conversion from a GW to 
co-propagating EM- waves is directly proportional to the background field energy density as well as the length of the 
interaction region, defined as the region occupied by the static magnetic field Bq. This conclusion holds as long as 
QED effects is negligible, and the length of the interaction region is smaller than the radius of curvature associated 
with the magnetic field energy density. Our aim here is to investigate to what extent the QED effects, associated with 
fields strengths approaching the Schwinger limit, modifies the energy conversion between GW:s and EM-waves. For 
this purpose we will still assume that the interaction region is smaller than the radius of curvature due to Bq, such 
that the interaction can be considered as taking place on a Minkowski background. 

As we will see, in addition to an EM-wave co-propagating with the monochromatic GW, with metric perturbation 
hx,+ = /ix,-i-exp [i{kz — cot)] and w = kc, a counter-propagating wave with the same frequency will also be induced. 
We thus make the ansatz B = JSoei -I- S'B{z) exp[—iujt] and E = (5E(z) exp [— zwi], where (5B and (5E includes both 
positive (along z) and negative propagating waves. Taking the curl of the pUj) and using ^ one obtains. 



e 



5B - V X (V X B) + ^oV X JQ = -^oV X JB + Moeojs, (13) 



to linear order, with the components of the polarization current Eq. (|lll) given by 

f = -l££B'Q^eoSE, , jI = -lE^Bl^d^m, , j| = 0. (14) 

^ IF cpo ^ 7F Mo ^ 



From Eq. ([12]) and Eqs. (IBip the gravitational contribution is found to be: 
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Bq dh+ .3 
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0, 



Jb = 0. 



(15) 
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Using Eqs. ([9]), ([T3|) . (|14p and (|15p we will next demonstrate that different EM wave polarizations couple to different 
GW polarizations. The result is most easily expressed in terms of the magnetic field components, and can then be 
written: 









SBi 
6B2 



k]^ Bdh^ exp [ikz] 



— + A:^ ) Sq/ix exp [ikz] , 



(16) 



where fc^ 



uj"^ / (c^ (1 + BqJff/if)) and k^ = w^ (l - B^jcg/if) /c^- As can be seen, all effects of the QED- 
vacuum polarization and magnetization is encoded in the effective wave-numbers fc^ and k^, that approach uj/c for 
cBq/Eci- -C 1. Note that Eq. (|16l) agrees with Ref. [I3, when the GW-coupling terms on the right hand sides 
are dropped [231 ■ The backreaction on the GW can be obtained by combining Eqs. ([6]) and (JA11|) . Whether or 
not this effect is important depends on the ratio of the excited wave energy density compared to the (pseudo) wave 
energy density of the GW. Roughly the scaling is as follows: For weak background magnetic fields (i. e. negligible 
QED effects), the excited wave energy density is limited by Wem '^ Bf/fiQ ~ (fcL)^|/i+,x pSo/^o , where k is the 
incident wave number and L is the length of the interaction region. As we will see in the next section, whenever QED 
effects are important, the excited wave energy is reduced compared to this scaling. Thus at most the the ratio of the 
excited wave energy to the GW (pseudo) wave energy density becomes Wcm/W^GW = i^(G/87rc^)_Bo/A'o- Whenever 
the interaction region is smaller than the background curvature due to the unperturbed magnetic field (as we have 
assumed above), this ratio is much smaller than unity, and hence the backreaction on the GW can be neglected. As 
a consequence, the approximation of "no GW back- reaction" will be employed in the next section. 



IV. A SPECIFIC EXAMPLE 



As a specific example we will now consider a boundary value problem, where the GW propagating in the 63- 
direction, enters the interaction region, given by —L/2 < z < L/2, which is the region where the external magnetic 
field B^ei is taken to be nonzero. The general solution to Eq. ([T6|) . for the interaction region —L/2 < z < L/2, is 



SB, 2 = Ti,2e'*=E 



i?i2e 



^ikt 



Cl 26 



ikz 



where Ci = fe^ BQh^/{k^ — fc^), C2 — {k^ + fe^)i?o^x/2(fc^ — k"^), and i?i^2 and Ti_2 are constants determined by 
the boundary conditions. This must be matched with the EM wave solutions with constant amplitudes outside the 
interaction region 

^5i,2 = /iV''" '^e (-^,-L/2), 



SBi,2 = fl 



Akz 



,z e (L/2, 00), 



at z — ±L/2. Furthermore, the electric fields must be matched as well. The relevant Maxwell equations are 




dzSBi 



dJBo 



Bo 



d^h 



zll'+ 



-^Cf^ft-x 



(17) 
(18) 



The matching of the electric field is done in the same way as that of the magnetic field to give four equations for four 
quantities, for each set of coupled polarizations. Solving these equations, the resulting amplitudes of the "refiected" 
and "transmitted" (or strictly speaking counter-propagating and co-propagating) EM- waves becomes 



J, _ Boh+ _,s (1 + 7;+) e"'+« + (1 - 77+) 6-'"+^ 

/l — — o — ^+^ 



2e" 



(1 -I- T]+) 6"*'^''+ - (r/+ - 1) e*'^i+ 



(19) 
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-irjj^O 
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(20) 



for the mode that couples to the plus-polarization. For the mode that couples to the cross-polarization we similarly 
obtain 



I^ = ^ {ril + 1) 



V 



^iSvx 



,-i&Vx 



(ryx - 1)^ e^^ix - (r?x + 1)^ e-*^''x 



and 
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(21) 



(22) 



Here we have introduced the notation 77+. x = ^b ^/^ and 9 = kL. An example of the magnetic profile (containing 
both the transmitted and reflected wave) is given in FiglT]for kL = 40 and cBo/Ea = 100. The expressions ([19])- 




FIG. 1: The wave profile for kL — 40 and cBo/Ecr = 100. The magnetized region lies between z = —10 and z = 10. 



(|22p contains all information about the energy conversion to the different EM-modes. However, to appreciate these 
results and the effects due to QED, we must first evaluate some results for the low-field limit when 77+^ x ^ 1- The 
squared coefficient |/i"|^, proportional to the energy density of the transmitted wave excited by the -I— polarization, 
then becomes 

o 1 2 

(23) 
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^1 -i 
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\ieL^ 


)site polarization. 
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\lk^L^ 



(24) 

Thus we see that the transmitted energy density is directly proportional to the background energy density. However, 
this behavior is dramatically changed when QED-effects are taken into account. The main reason is that the EM 
wave dispersion relation is changed in the interaction region (that makes 774-. x deviate from unity) which in turn 
detunes the excited wave with the GW. The consequence for the transmitted wave excited by the x -polarization 
is depicted in Fig. [21 for kL = 20 and kL = 100. The steady increase in the absence of QED is replaced by an 
oscillatory behavior, mainly due to the detuning of the GW and EM wave dispersion relation. Note that we here 
have normalized the transmission coefficient with \h+ ^^\^ B^k'^ L'^ , such that the coefficient without QED-effects is 
represented by a straight line. For a longer interaction region, a smaller mismatch of dispersion relations are needed 
for the phase difference to accumulate, and hence the curve with the lower value of kL (kL = 20) needs a much higher 
field strength before significant QED-effects are seen. A similar point is illustrated by Fig. [3] that depicts the energy 



density for the co-propagating mode excited by the +- polarization. Note that the energy conversion to this EM-mode 
is much less affected by the QED effects. The reason is that the QED-modification of the EM dispersion relation 
effectively saturates at a value cBq/Ecy ~ 10. Accordingly we have chosen higher values of kL, namely kL = 2000 and 
kL = 20000, which is needed in order to see the deviation from the classical behavior induced by QED. In addition to 
the co-propagating EM modes there are also counter-propagating EM- waves. From a practical point of view, these are 
much less significant, since the counter-propagating modes are always non-resonant with the source GW, and hence 
the energy density of these modes does not systematically increase with a larger interaction region, i.e. increasing 
kL. From a more theoretical point of view, an interesting effect can be seen in the coefficients ((T^ and (l?T|) . however. 
Without QED-effects, the -F-polarization does not cause a back-scattered wave, independent of the value of kL, as 
seen by ([T5|) when letting rj^ —^ 1. However, the situation for the x -polarization is different, as we find a finite but 
small counter-propagating mode from (|2ip also in the limit rj^ ^ ^■ 
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FIG. 2: Normalized energy density of the copropagating EM-wave, excited by a cross-polarized GW, as a function of background 
magnetic field strength for 9 — 20 (dashed line) and 6 — 100 (solid line) compared to the non-QED-case (dotted line). 
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FIG. 3: Normalized energy density of the copropagating EM-wave, excited by a plus-polarized GW, as a function of background 
magnetic field strength for 9 — 2000 (dashed line) and 9 = 20000 (solid line) compared to the non-QED-case (dotted line). 



SUMMARY AND CONCLUSION 



In this paper we have studied the interaction between GW:s and EM- waves in the presence of a strong static magnetic 
field Bq, using the Heisenberg-Euler lagrangian in order to take QED vacuum polarization and magnetization into 
account. The high-frequency approximation has been applied to zeroth order, i.e. all effects of the background 
curvature has been neglected, which is permissible if the spatial extension of the interaction region is much smaller 
than the radius of curvature. The specific boundary conditions considered is an incoming GW incident on a static 



magnetic field with a given extent L, which give raise to an excited EM-wave in the same direction as the GW, as well 
as one propagating in the opposite direction. The role of the QED effects is twofold: Firstly, the coupling strength 
between the GW:s and the electromagnetic waves are modified (as described by the coefficients of the right hand side 
in Eq. ([T5| '). Secondly, the change in phase velocity (< c) of the EM- waves induced by the vacuum polarization, as 
described by the expressions k^ and kj,, destroys the perfect resonance with the gravitational source wave, which 
gives a saturation of the possible energy conservation at a finite value of L. These effects are similar in principle 
for the hx- and /i+-polarizations (which couples to different EM-polarizations), and the dimensionless parameter 
{cBo/EcY)'^kL need to reach [cBq/ Eci)^kL ^ 10^ in order for QED effects to be important in both cases. However, 
since the QED-modification of the EM-mode excited by the /i+-polarization saturates at a value cBq/Eci ^ 10, a 
much higher value of kL is needed for the QED-effects to be significant in this case. 

The problem considered here has been highly idealized and has mainly been motivated by a theoretical interest 
to study GW and EM-wave interaction in a strong field environment, allowing for field strengths larger than the 
Schwinger critical field E^- However, we would like to point out that there is a certain astrophysical relevance of the 
problem, as the effect of QED-detuning is found to be of significance for field strengths Bo — 3-Ecr/c ~ 10^°T (see e.g. 
Fig [3]), a value that has been observed at magnetar surfaces [2l|, although a high GW frequency would be required. 
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Appendix A: Strong field vacuum polarization and magnetization parameters 

With only a strong magnetic field present the quantities 7j^,7g,7j^j^,7ge and jjrg can be determined analytically, 
see Ref. pjj. The resulting expressions for these QED-parameters are 



7e = 0, 
Ij^g = 0, 
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/io Stt^o 
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where a = . ^ > is the fine structure constant, h ~ j^, T{h) the gamma function, ij^^h) the digamma function and 
C'(— l,/i) the first derivative of the Hurwitz zeta function with respect to its first argument. 

Furthermore, in the absence of a strong electric field we can calculate the integral in the Lagrangian ^ analytically. 
Since there is only a strong magnetic field present we have 6 = 0. Thus, to compute the integral in Eq. (IT|), we expand 
the integrand and take the limit as 6 — >■ 0, thereby obtaining 



"■'' -eEgTs/c 



(cdS) 

(eas) coth (eas) 1 
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By changing the variables such that eas — z, dividing the integral into three parts, altering the integration path and 
using the regulator z*^ we obtain 



/ = (eay 
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Since Ecr/ca — 2h we find the first, second and third part of the integral to be 
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dz e-^'^'-z^--^ = - [ + /i^ - 2/i^ In /i + (1 - C - hi 2)2h^ 



respectively, where C is Eulers constant. With Eas. (jA4p . (|A5P and (jA6[) we can now rewrite Eq. (jA2l) as 

/ = {eaf I i [1 - In 2 - In (2/i)] + /i^p In /i - 1] - 2/i In /i - 4C'(-1, /i) j , 

and thus the Lagrangian ^ becomes 

C = — ^-F- -^^|i[l-ln2-ln(2/i)] + /i2[21n/i-l] -2h\nh - AC'{-l,h)\ - A^j" . 
^0 27r/xo [3 J 

Since we have only a magnetic field, G — holds, and the energy-momentum tensor associated with the La; 
(IA8|) becomes 

Next we proceed by expanding the energy- momentum tensor (jA9p . The first order contribution becomes 

where 
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and ho ~ Ecr/2cBo, so the relevant energy-momentum tensor terms in Eq. ([6]) becomes 

STii - ST22 = BqSjf -2jFBQdBi, 
ST12 — ■yFBoSB2- 



(All) 



Appendix B: Ricci-rotation coefficients 

The Ricci-rotation coefficients of a Minkowski spacetime perturbed by a GW propagating in the ea-direction 
expressed in the tetrad (O is given by 

7?i = -722 = 7oi = -7o2 = ^^+' 



1 ■ 
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to first order in h 



+,x- 



[1] J. Moortgat and J. Kuijpers, A&A 402, 905 (2003). 



[2] M. Marklund, G. Brodin and P. K. S. Dunsby, Astrophys. J. 536, 875 (2000). 

[3] H. Isliker, I. Sandberg, L. Vlahos, Phys. Rev.' D 74, 104009 (2006). 

[4] D. Papadopoulos, N. Stergioulas, L. Vlahos and J. Kuijpers, A&A 377, 701 (2001). 

[5] A. Kallberg, G. Brodin and M. Bradley, Phys. Rev. D 70, 044014 (2004). 

[6] Yu G. Ignat'ev, Phys. Lett. A 320, 171 (1997). 

[7] G. Brodin and M. Marklund, Phys. Rev. Lett. 82, 3012 (1999). 

[8] G. Brodin, M. Marklund and P. K. S. Dunsby, Phys. Rev. D 62, 104008 (2000). 

[9] G. Brodin, M. Marklund and M. Servin, Phys. Rev. D 63, 124003 (2001). 
[10] A. B. Balakin, V. R. Kurbanova and W. Zimdahl, J. Math. Phys., 44, 5120 (2003) 
[11] M. Servin and G. Brodin, Phys. Rev. D 68, 044017 (2003). 

[12] M. Servin, G. Brodin, M. Bradley and M. Marklund, Phys. Rev E, 62, 8493 (2000). 
[13] D. Papadopoulos, Class Quantum Grav. 19, 2939 (2002). 

[14] M. Marklund, P.K.S. Dunsby, and G. Brodin, Phys. Rev. D 62, 101501 (2000). 
[15] H. J. M. Cuesta, Phys. Rev. D 65, 64009 (2002). 
[16] J. Moortgat and J. Kuijpers, MNRAS 368 , 1110 (2006). 
[17] M. Marklund and P. K. Shukia, Rev. Mod. Phys. 78, 591 (2006). 

[18] G. Brodin, L. Stenflo, D Anderson, M. Lisak, M. Marklund and P. Johannisson, Phys. Lett A, 306, 206 (2003). 
[19] J. Lundin, Europhys. Lett., 87, 31001 (2009). 

[20] S. R. Valluri, D. R. Lamm and W. J. Mielniczuk, Can. J. Phys., 71, 389 (1993). 
[21] C. Kouveliotou et al, Nature 393, 235 (1998). 
[22] J. Schwinger, Phys. Rev. 82, 664 (1951). 

[23] W. Dittrich and H. Gies, Probing the Quantum Vacuum (Springer- Verlag, Berlin) 2000. 
[24] G. F. R. Ellis and H. van Elst, Cosmologtcal m,odels, Theoretical and Observational Cosm,ology, ed. M Lachieze-Rey 

(Dordrecht: Kluwer) (1999). 
[25] It should be noted that although the interaction Eq. (|16|) can be expressed solely in terms of 5B\ and 5B2 , the EM- wave 
polarization in the presence of strong QED-effects is nontrivial. For a more complete description, see e.g. Ref. [IS] . 



